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Abstract. We study the biased diffusion of particles moving in one direction under the action of a constant 
force in the presence of a piecewise linear random potential. Using the overdamped equation of motion, we 
represent the first and second moments of the particle position as inverse Laplace transforms. By applying 
to these transforms the ordinary and the modified Tauberian theorem, we determine the short- and long- 
time behavior of the mean-square displacement of particles. Our results show that while at short times the 
biased diffusion is always ballistic, at long times it can be either normal or anomalous. We formulate the 
conditions for normal and anomalous behavior and derive the laws of biased diffusion in both these cases. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.10.Gg 
Stochastic analysis methods (Fokker-Planck, Langevin, etc.) - 02.50.-r Probability theory, stochastic pro- 
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1 Introduction 

The biased diffusion, i.e., diffusion accompanied by the 
directional transport of diffusing objects (particles), ex- 
ists in many systems subjected to external force fields. In 
systems with quenched disorder it can be naturally de- 
, scribed by the Langevin equation in which disorder is ac- 
■ counted by time-independent random potentials and ther- 
mal fluctuations are modeled by white noise. Because of its 
relative simplicity and efficiency, the Langevin-based ap- 
proach provides an important tool for studying the trans- 
port properties in disordered media [T]. In particular, wit- 
hin this framework a number of effects arising from the 
joint action of quenched disorder and thermal fluctuations, 
including some features of biased diffusion, has been suc- 
cessfully studied for particles moving under a constant 
force in a one-dimensional random potential [2,3,4,5,6,7, 

■US]. 

Since in the most cases considered earlier the distri- 
bution of the random force which corresponds to a given 
random potential has unbounded support, particles can- 
not be transported to an arbitrary large distance if ther- 
mal fluctuations are absent. Put differently, in the noise- 
less case particles remain localized in a finite region at 
all times. It is clear that derealization can occur only if 
the above mentioned distribution has bounded support. 
At this condition, the biased diffusion can be caused by 
a time-periodic external force or a constant one. In both 
these cases the diffusive behavior of particles results solely 
from quenched disorder, but the directional transport has 
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different nature. Indeed, a time-periodic force induces the 
directional transport due to the ratchet effect which exists 
in random potentials of a special class, i.e., ratchet poten- 
tials with quenched disorder [l0l[TTl[T2l[T3] . In contrast, a 
constant force (if it exceeds a critical value) induces the 
directional transport in arbitrary potentials due to the di- 
rect action on particles. In this case particles move only 
in one direction and, as a consequence, the completely 
anisotropic case of biased diffusion, when the probabil- 
ity of motion along and against the external force equals 
1 and 0, respectively, is realized. Some features of this 
diffusion exhibiting normal behavior were considered in 
P^lfTSlfTB] within a general approach based on the equa- 
tion of motion of diffusing particles. At the same time, 
anomalous regimes of this diffusion were studied only in 
the framework of continuous time random walk [17] , which 
adequately describes the long-time case. 

In this paper we develop a unified approach for the 
study of biased diffusion of particles moving under a con- 
stant force in a piecewise linear random potential. It is 
based on the overdamped equation of motion and gives a 
possibility to get the diffusion laws for both short and long 
times. The paper is organized as follows. In Section [5] we 
describe the model and derive the probability density of 
the particle position in terms of the probability density of 
the residence time. In Section [3] we represent the first two 
moments of the particle distribution as inverse Laplace 
transforms. The biased diffusion at short and long times 
is studied in Sections U and [SJ respectively, using the ordi- 
nary and the modified Tauberian theorem for the Laplace 
transform. Specifically, we obtain the law of diffusion at 
short times (Sections 2]), formulate the conditions for nor- 
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Fig. 1. Sample path of the piecewise constant random force 
g(x) that corresponds to a given realization of the piecewise 
linear random potential (UJ). 



mal and anomalous diffusion at long times (Section I5.1[) . 
calculate the diffusion coefficient and analyze its depen- 
dence on the driving force (Section 15.21) . and derive the 
laws of anomalous diffusion (Section 15.31) . Finally, in Sec- 
tion [6] we summarize our results. 



2 Probability density of the particle position 

We consider the unidirectional motion of a particle which 
occurs under the action of a constant driving force /(> 0) 
in a piecewise linear random potential U(x). At x > we 
define this potential as follows: 

U{x) = -{x-nl)g (n) + U(nl), x € [nl,nl + I), (I) 

where n = 0,1,... and the slopes are assumed to 
be independent random variables having the same proba- 
bility density function u(g). Next we suppose that this 
density function is symmetric, i.e., u(—g) = u(g), and 
has bounded support, i.e., g G [— go, go]- If the condition 

U(nl) = -ZEm=U (m) holds for n > 1 ( U (°) is an arbi- 
trary constant which can be chosen to be zero) then the 
potential U(x) is continuous and the corresponding piece- 
wise constant random force, g(x) = —dU (x)/dx, takes the 
form g(x) — if x belongs to the nth interval [nl, nl + l) 
(see figure [T]). Thus, neglecting the inertial effects, we can 
describe the particle dynamics by the overdamped equa- 
tion of motion 

vX t = f + g(X t ), (2) 

where X t (Xq = 0) is the particle position and v is the 
damping coefficient. 

According to the above assumptions, at / G (0,<?o) 
particles cannot be transported to an arbitrary large dis- 
tance in the positive direction of the axis x. On the con- 
trary, they are stopped at some distance L — IN from the 
origin, where N is the number of that interval which is 
characterized by the conditions / > — g^ (for all n < N) 
and / < —g^ N ' This number depends on the sample paths 
of g(x) and its mean (N) (the angular brackets denote 
an average over these sample paths) can be easily calcu- 
lated by introducing the probability / = J_ dgu(g) that 

g(x) < /. Indeed, since the random forces g^ are statis- 
tically independent on different intervals, the probability 



W n that N = n can be written as W n = I n f?° dgu(g), 
i.e., W n = I n — I n+1 . With this expression for W n we 
obtain (N) = YlnLi n W n = and the use of 

the geometric series formula E^Lo^" = — -0 yields 
(N) = 1/(1 — I). As a consequence, the average distance 
(L) to the point where particles are stopped is given by 



(L) = l 



1 -I 



(3) 



Thus, if / < go, i.e., / < 1, then the final state of particles 
is their localization at a finite average distance ([3]). We 
note also that if the probability density u (g) has no 5 
singularities at g = ±go, i.e., u(g) does not contain the 
terms proportional to 5(g — go) and 5(g + go), where S(-) 
is the Dirac S function, then I — > 1 and (L) — > oo as 

/ -> Sta- 
in contrast, if / > go then f+g(X t ) > and, according 
to the motion equation @, X t — > oo as t — » oo. As it 
follows from ((3J), the same behavior holds for / = go as 
well. In other words, at / > go particles can be transported 
to an arbitrary large distance along the positive direction 
of the x-axis. In this cases the distribution of particles on 
the positive semi-axis x is described by the probability 
density 

P(x,t) = (6(x - X t )) (4) 

that Xt = x. Assuming that X t — nl + x[ n ^ and x\ G 
[0,1), it is convenient to rewrite the probability density of 
the particle position in the form 



P{x,t) = Y j P n {x,t) 



(5) 



with P n (x,t) = (S(x — nl — Xt)). In order to calculate 
P n (x,t), let us first introduce the residence time t'"J of 
a particle in the nth interval [nl,nl + I), i.e., time that 
a particle spends moving from the point nl to the point 
nl + l. Since in the overdamped regime the motion occurs 



with a constant velocity = l/T (n '> , we obtain x[ n ^ — 
Kt - ElTio tW)/tW if t ~ £r=o rM G [0, rW) and 



Pn(x,t) 



nl 




(G) 



According to (0), the residence time is given by t^ 1 ' = 
vl l '(/ ' + gW) . Because of the properties of g^ , these times 
in different intervals are statistically independent and dis- 
tributed with the same probability density p(r). Taking 
into account that r G [r m i n ,r max ] with 



vl 



vl 



f + go' 



f - go 



(7) 



from the condition p(r)|dr| = u(g)\dg\ we can express the 
probability density p(r) of the residence time through the 
probability density u(g) of the random force: 



p{r) 




min? ''"max J 

otherwise. 



(8) 



S.I. Denisov, E.S. Denisova, H. Kantz: Biased diffusion in a piecewise linear random potential 



3 



In turn, the quantities P n (x,t) can be expressed thro- (Res > 0) to both sides of this formula. Using the convo- 

ugh p(r). Specifically, from the definition ([6]) at n — we lution theorem for the Laplace transform, £{v(t)*w(t)} = 

obtain v s w s , this yields 

P (x,t) 



(16) 



It is also not difficult to see that P n {x,t) at n > 1 is given 
by By changing the order of integration in the Laplace 



transform of the integral term, we obtain 
Jo Jo c\ ( dT(n + t/T)p(T)\ = [ dTp(r) 



P n (x,t) = / dT p{T ) / dnpirx) 



xdr n _ip(T„_i) / dT n p(T n ) f r , 

-E:=o t » x / dt(n + t/r)e- st = n i— ™ ~^+F s , (17) 



■s 



X<5 



\ m— 



(io) 



where 

1 f°° 1 - e~ ST 
F 5 = -/ p(r). (18) 



The last two formulas can be written in a more compact 

form using the notation v(t) * w(t) for the convolution of Therefore, taking into account that E^Lo^™ = V(l ~~Ps) 
two functions v(t) and w(i), and E^Li n Ps = Pa/ (1 — Ps) 2 , formula (|16|) reduces to 

v(t)*U,(t)= f dTv(T)w(t-T). (11) ( Xt ) s =Z-^. (19) 

JO 1 - Ps 

With this definition, we can introduce the n-fold convolu- applying to © the inverse Laplace transform 

tion of the probability density of the residence time as fol- defined as 
lows: p* n (t) = p*™- 1 ^) *p(t). Assuming also that p*°(t) = 

5(t), we obtain P (x,t) = p*°(t) * P Q (x,t), P n (x,t) = _ i r c +^ 

p* n (t) * P (x - nl,t), and so the probability density of h(t) = C {h s } = — / dse st h s (20) 

the particle position can be represented in the form 

oo (the real number c is chosen to be larger than the real 

P(x,t) = Y^p* n (i) * Po(x — nl,t). (12) parts of all singularities of h s ), we find the first moment 

n= o in the form 

We note that there is a possibility of further simplifica- (Xt) =IC x < - > ■ (21) 

tion of P(x, t) by using the Laplace transform method. ^ Ps ' 

However, the representation (|12p is quite sufficient for our The second moment, 
purpose, i.e., finding the moments of X t . 



^ i>oo 

(X?) =l 2 J2 P* n (t) * / dr(n + t/rfp(r), 

„-n Jt 



( 22 ) 

3 Moments of the particle position 

The fcth moment of the particle position is defined in the f n a1 ^ be determined by the method of Laplace trans- 

i * form. Ihe similar calculations lead to 

usual way: 



J™dxx k P(x,t) (k = l,2,...). (13) (X*) a = l 2 f^Cy™fr(n + t/T) 2 p(T)y (23) 



Substituting (fT2l) into (fT3]) and using ([9]) , for the first mo- wne re 
mcnt we obtain 



°° poo 

{X t ) = lJ2p m (t) * / Mn + t/r)p(T). (14) 

Because of the convolution structure of the summands in 
HH), it is reasonable to apply the Laplace transform de- 
fined as with 



Cl I dr(n + t/Tfp{r) \ = n 2 ?—El + 2n ^ Ps 



P - + G s (24) 



2 r,l-(l + sr)e- 



h s =C{h{t)}= dte' st h(t) (15) G s = ^ dr ' p(r). (25) 
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Carrying out summation over n in equation (|23[) [for this the first two moments of Y t can be written as 
purpose we use the above expressions for the series and 



the formula Y^=i n2 P" = Ps(^ +Ps)/(1 -Ps) 3 }, we obtain /y t ) = l£~ 



and so 



(X?) 



^PsF s 



(l-Ps)G s 



{i-PsY 



(26) 



S{l-Ps) 



, (y 2 ) = i 2 £- 



Ps+Ps 

s(l- Ps ) 2 



The difference in the time dependence of the variances 
a 2 x {t) and c4(i) = (Y t 2 ) - (Y t ) 2 will be clarified in the 
next sections. 



(X, 



2 Ps F s + (l- Ps )G s 
(l-Ps) 2 



(27) 



As it will be shown below, the representations (|2"Tj) and 
(|27|) of the moments (Xt) and (X 2 ) are very useful for 
finding the short- and long-time behavior of the mean- 
square displacement of a particle, or the variance of the 
particle position, defined as 



(x 2 ) - (x t f 



(28) 



In conclusion of this section, we briefly discuss the dif- 
ferences between our model and the continuous-time ran- 
dom walk (CTRW) models. Introduced many years ago 
[18] , the CTRW approach has become one of the most im- 
portant tools for the study of anomalous transport (see 
e.g. [I], [1*12011211125 1125] ). In the simplest version of the 
CTRW is assumed that the random waiting time, i.e., time 
that a walking particle waits before it jumps to another 
position, and the random jump size are statistically in- 
dependent. Since equation ([2]) describes the continuous 
motion of a particle, this so-called jump model is in gen- 
eral inapplicable to our case. Due to physical reasons, 
much attention has also been paid to the case in which 
the waiting time (considered as the motion time) and the 
jump size (considered as the passed distance) are coupled 
in such a way that a walking particle moves with a con- 
stant velocity 24 , 25 ,26,27, 28 . However, because of the 
last condition, this velocity model also does not provide 
an adequate description of the biased diffusion in a piece- 
wise linear random potential. From this point of view, 
the CTRW model [29] in which a walking particle moves 
with a random velocity could be the most appropriate. 
But in [20J is assumed that the random velocity and mo- 
tion time are statistically independent, while in our case 
the particle velocity »W = (/ +g^ n ')/u and the residence 

time t^") = vl/{f + g^ n ') are perfectly dependent, i.e., 

v (n) T (n) = i 

Thus, none of the CTRW models completely describes 
the considered type of biased diffusion. Nevertheless, if 
the waiting time is associated with the residence time and 
the jump size with I then the long-time behavior of X t 
can be adequately described by the CTRW jump model 
[T7] . It is therefore reasonable to compare the behavior 
of X t with that following from the unidirectional CTRW 
on a semi- infinite chain of period I. Within this approach, 
the particle position is described by a discrete variable 
Y t — lN(t), where N(t) is the random number of jumps 
up to time t, and the waiting time is characterized by the 
same probability density p(r). Using the moments of N(t) 
that are known from the theory of CTRWs (see e.g. [H]), 



4 Short-time behavior of the variance 

The calculation of the inverse Laplace transforms in for- 
mulas (|2"T]) and (|27D is a difficult technical problem. For- 
tunately, in the short- and long-time limits this prob- 
lem can be avoided. Such a possibility provides the cel- 
ebrated Tauberian theorem for the Laplace transform [19, 
ISO] which states that if h(t) is a monotonic function and 



h, = C\h(l)} [ 1)1 



as s — > (s — > oo) then 

h{t)=C- 1 {h s } 



r(p) 



L(t) t p - 



(30) 



(31) 



as t — > oo (t —> 0). Here, p G (0, oo), r(x) is the gamma 
function, and L(t) is a slowly varying function at infinity 
(zero), i.e., L(Xt) ~ L(t) as t — !• oo (t — > 0) for all A > 
0. It should be stressed that, in contrast to the Laplace 
transform (TT5]) . the parameter s in (j3"0"]) is assumed to be 
a positive real number. We note also that if s — >• then 
the Tauberian theorem is valid in a more general case, 
namely, if h{t) is ultimately monotone, i.e., monotone on 
the interval (£', oo) for some t' > 0. 

Thus, for finding the short-time behavior of the mo- 
ments (Xt) and (X 2 ) at / > go we should determine 



the leading terms of the Laplace transforms (Xt) 

^2 



and 



(X 2 ) s as s tends to infinity. Taking into account that 
t <= bmin, Tmax] and p s — > as s — > oo, from (fT9]) and 
(1261) in this limit we obtain 



I 



dr^ (X 2 ) t 



2P 

B 3 



dr 



P( T ) 
t 2 ' 
(32) 

Using ((HJ and the symmetry property of the probability 
density of the random force, u(—g) = u(g), the integrals 
in can be expressed as 



dr 



P(r) _ f 

T vl ' 



(It 



p(r) _P + °l 



(ulf 



(33) 



where a 2 = fH° go dgg 2 u(g) is the variance of the random 
force g(x). Therefore, as follows from the Tauberian the- 
orem, the short-time behavior of the first two moments is 
described by the asymptotic formulas 



(X t )~ J -t, (Xt) 
v 



(34) 
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They show that at small times (t <C vlj f) and f > go the with positive parameters a and a. According to it t 2 < oo 



biased diffusion is always ballistic 
tional to t 2 : 



i.e., &x(t) is propor- 



2 L ■ 



(35) 



We note that the diffusion law (l3"5j) also follows from the 
motion equation ©. Indeed, the first two moments of the 
particle position X t ~ t(f + g^°')/v, which is the solution 
of this equation at t — > 0, are the same as in ([Ml) . 

For finding the short-time behavior of the moments 
of the discrete variable Y t the Tauberian theorem cannot 
be used. The reason is that (Y t k ) s at s — > oo tends to 
zero more rapidly than any positive power of 1/s (since 
p s oc e~ STmin as s -> oo). Nevertheless, exact conclu- 
sions about the short-time behavior of (Y t ) can be drown 
directly from the definition (Y t k ) = l k Y^=i n k V(n, t), 
where V(n,t) is the probability that N(t) = n. Indeed, 
accounting for the fact that V(n,t) — p(t) * V(n — l,t) 
(n > 1) and V(0,t) = f°° drp(r) [T5], we make sure (see 
also © and (JTlJ) ) that if t < r min then 7>(0,t) = 1 and 
V(n,t) = for all n > 1. Hence, at t < r m - m all moments 
of Y t are equal to zero and so the short-time behavior of 
Y t is qualitatively different from that for X t . 



5 Long-time behavior of the variance 

5.1 Conditions of normal and anomalous diffusion 

According to ([5]), the mth moment of the residence time, 
r m = /„ dTT m p(r) (to = 1,2,.. .), can be written in the 
form 



(yi) 7 



^ d9 (f+ 9 y 



(36) 



Since the probability density u(g) is properly normalized 
on the interval [— go, go], i-e., f_° o dgu(g) = 1, the above 

formula leads to the condition r™ < {vl) m /{f - g ) m . It 
shows that if / > go then all moments r m are finite. Since 
in this case r 2 < oo, the classical central limit theorem 
for sums of a random number of random variables |31j is 
applied to A t and suggests that cr\{t) oc t as t — > oo. In 
other words, at / > go the biased diffusion of particles is 
expected to be always normal. 

In contrast, if / = go then the residence time varies 
from T m ; n = vl/2go to r max = oo, and t 2 can be either 
finite or infinite depending on the asymptotic behavior 
of p(r) as r — > oo (i.e., asymptotic behavior of u(g) as 
g — > —go)- While in the former case the biased diffusion is 
also expected to be normal, in the latter (when r 2 = oo) 
the above mentioned central limit theorem becomes inap- 
plicable and one may expect that the biased diffusion is 
anomalous. Next we assume that at / = go the probability 
density p(r) of the residence time is characterized by the 
asymptotic formula 



if a > 2 and t 2 = oo if a € (0, 2]. The parameters a and a 
depend on the explicit form of the probability density u(g) 
of the piecewise constant random force g(x) and are in 
general not independent. In particular, for the symmetric 
beta probability density 



r(P + i/2) 
goV^r(P) 



So 



0-i 



(f3 > 0) we have a = (3 and 

2vl 
9o 



r(a + 1/2) 
2^r{a) ' 



(38) 



(39) 



Thus, the biased diffusion is expected to be normal 
if / > go or / = go and a > 2, and anomalous if / = 
go and a € (0,2]. In the next two sections we validate 
these criteria and derive the laws of normal and anomalous 
diffusion. 



5.2 Normal biased diffusion 



5.2.1 Modified Tauberian theorem. 

The ordinary Tauberian theorem formulated in (|30[) and 
(f3"Tj) permits us to find only the leading terms of the 
asymptotic expansion of h(t) at t — > and t — > oo. There- 
fore, in this form it can be used for determining the long- 
time behavior of the variance a\ (t) only in the case when 
the leading terms of (X t ) 2 and (X?) are different (see sec- 
tion But in order to find u\{t) in the case with 
{X t } 2 ~ (A 2 ), we need to know at least two leading terms 
of each of the asymptotic expansions of (X t ) and (A 2 ). 
With one exception (see section [5. 3. 3|) . these terms can be 
found by applying the modified version of the Tauberian 
theorem. This version follows directly from the ordinary 
one by replacing h s by h s — r j s r ' and using the exact re- 
sult L^il/s^} = ^- 1 /f(r?) [32]. According to this, we 
formulate the modified Tauberian theorem as follows: If 
h(t) — rf 1 ^ 1 / T{rj) is ultimately monotone and 



\s sp 



(?/ > p > 0) as s — > then 
hit) 



r( v ) 



r{p) 



L(t) tP- 1 



(40) 



(41) 



p(r) 



(r 



oo 



(37) 



as t — > oo. 

It is important to emphasize that the terms rt v ~ l / r(rj) 
and L(t)t p_1 / T(p) represent the first two terms of the 
long-time expansion of h(t) only if rj > p. The reason is 
that in the opposite case (when rj < p) the second term 
of the asymptotic expansion of £~ 1 {L(l/s)/s p } at t — > oo 
may not be negligible in comparison with C~ 1 {r/s v }. In 
order to illustrate this fact, we assume that L(l/s) — 
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ln(l/s) and consider the Laplace transform of the follow- 
ing form: 

^ + (42) 

Since £- 1 {ln(l/s)/s''} = tP- l []nt-ip(p)]/r(p) [32], where 
ip(p) is the digamma function, i.e., the logarithmic deriva- 
tive of the gamma function r(p) , the inverse Laplace trans- 
form of h s yields an exact result 



h(t) 



t'i- 



r(v) r(p) 



pnt-v(p)]. 



(43) 



Keeping in the long-time limit two leading terms of h(t), 
we obtain 

( t^/rirj) + tp- 1 In t/r(p), V > p 

h(t) ~ \ t"- 1 [lnt-V(?7) + l]/r(77), n = p (44) 

[tp- l [hit-^(p)}/r(p), V<P . 

By comparing the asymptotic formulas (|4"Tj) and ([4"4l we 
make sure that at r/ > p and t — > oo two leading terms of 
h(t) can indeed be determined from the modified Taube- 
rian theorem. In contrast, if 77 < p then for finding two 
leading terms of the long-time expansion of h(t) it is nec- 
essary to go beyond the Tauberian theorem. 

5.2.2 Properties of the diffusion coefficient. 

We define the coefficient of biased diffusion in the same 
way as for ordinary diffusion: 



D h = lim -Z- 

t—>oc Zt 



(t) 



(45) 



As it follows from the modified Tauberian theorem and 
the definition ([2"5]) , for determining the law of biased dif- 
fusion, i.e., the long-time behavior of the variance cr^-(i), 
we should find two leading terms of each Laplace trans- 
form, (X t ) s and (X 2 ) s , in the limit s — > 0. Assuming 
that the first two moments of the residence time, r and 
r 2 , exist, the Laplace transforms p s and F s at s — > 
can be represented by the following asymptotic formulas: 
p s ~ 1 — sr + s 2 t 2 /2 and F s ~ 1/s — r/2. Substituting 
them into (11911 . we obtain the desired result 



I 1 



T 1 



2f l 



(46) 



Similarly, taking also into account that G s ~ 1/s as s — > 0, 
from (|26[) we find with the required accuracy 



/V2\ 212 1 o,2 t2 - ?2 1 



(47) 



The asymptotic formulas (|46l) and (|4T| are particular cases 
of (l40l) in which L(l/s) is a constant. Therefore, the mod- 
ified Tauberian theorem leads to the following expressions 
for the moments (X t ) and (Xf) in the long-time limit: 



(Xt. 



T 2_ T 2 



I 

-t + l——^ 

t 2t 2 



l l 2 

rj i I 2/ 1 1 



(48) 



which in turn determine the law of biased diffusion 

_ Z2 _ ^2 



(TatC*) ~ ^ — =5 — *• 



(49) 



Thus, it follows from (1451) and (|49[) that the coefficient 
of biased diffusion can be written as 



Dh = l 2 - 



2r 3 ' 



(50) 



Using ([36"]) and the Jensen inequality [30] , one can see that 
for all probability densities u(g) the condition Db € [0, 00] 
holds. However, if / > go then the moments r and t 2 are 
finite (see section IBTTj) and nonzero (see formula ([36])), i.e., 
Db 7^ 00 ■ On the other hand, representing r 2 as 

r 2 = [ -J- I I dgdg'u(g)u(g') 



and t as 



r 2 = (vlf 



(f + g) 2 (f + g'f 



90 r d gdgi tt <*W> 

So So 



(/+. 9) (/ + .?') 



(51) 



(52) 



we obtain 



r^-r 2 



(vlf 



2 /■So /-So 



dgdg'u(g)u{g') 



So •'—So 

A2 



(/ + .g) 2 (/ + .g') 2 



(53) 



The last expression shows that always r 2 — r 2 > 0, and 
the condition t 2 = t 2 (D^ = 0) holds only if u(g) = 5(g) 
(we took into account the condition u(— g) = u(g)). But 
the probability density u(g) = 5(g) corresponds to the 
trivial case when the random force g(x) is absent. Ex- 
cluding it from consideration, we obtain 7^ 0. Thus, 
if / > go then D\ } g (0, 00), i.e., the biased diffusion is 
normal: cr\{t) ~ 2D^t. It is clear from the above analysis 
that at / = go and a > 2 the biased diffusion is also nor- 
mal. In contrast, if Z?b = or — 00 (these conditions 
are realized only if / = go and a € (0, 2], see section I5~31) 
then the biased diffusion becomes anomalous: It is slower 
than normal in the former case and is faster in the latter 
one. We note, however, that the conditions — and 
Z?b = 00 only indicate the existence of anomalous behav- 
ior and cannot be used for finding the laws of anomalous 
diffusion. 

Next we consider the dependence of D\> on the external 
force / and statistical characteristics of the random force 
g(x). Let us first determine Z?b in the case of large external 
force when / 3> go- By expanding the integrand of (I36[) in 
powers of g/f and integrating over g with the condition 
u{— g) = u(g), we obtain 




V 2 



(54) 
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Fig. 2. The reduced coefficient of unidirectional diffusion 
D — D^v/lgo as a function of the dimensionless driving force 
q — f/go- The monotonic dependence of D on q is determined 
from formula (|58|l which corresponds to the uniform proba- 
bility density u(g), and the non-monotonic one follows from 
formula ()63|) which corresponds to the probability density (|61l) . 



as / — > oo. The substitution of these asymptotic expres- 
sions into ([501) leads to the universal, inversely propor- 
tional dependence of Dh on / for large values of /: 



2v /■ 



(55) 



If / > go and the probability density u(g) is given by 
then, according to (1361) , the moments of the residence 
time can be represented in the form 



vl 



f - go 



F[ m,(3;2(3;- 



2go 
f -go 



(56) 



where F(m,/3;2/3;z) = 2F\(m, f); 2/3; z) is the hyperge- 
ometric function. In some cases the first two moments, 
t and r 2 , can be expressed in elementary functions. In 
particular, it is possible for j3 — 1, i.e., if u(g) corre- 
sponds to a uniform distribution of the random force g(x) 
on the interval [— go, go]- In this case, using the condi- 
tions [33] F(l, 1; 2; z) = -z" 1 ln(l - z), F{m, (3; 2/3; z) = 
F(/3, m; 2/3; z) and F(l, 2; 2; z) = (1 - z)' 1 , we get 



vl , q + 1 
In , 

250 q - 1 



go 



and, as a consequence, 



D h = 



Igo 



q 2 - 1 



-In 



2 9 



In 



-3 9- 



9-1' 



(57) 



(58) 



where q — f/go > 1 is the dimensionless external force. 
The coefficient of biased diffusion (1551) is a monotonically 
decreasing function of q (see figure [2]) which tends to in- 
finity as^l,A)~ (2lg /v)(q- l)" 1 ! hr 3 (<7 - 1)|, and 
approaches zero as q —> 00, D\, ~ (lgo/6v)q . Since for 
the uniform distribution o 2 g = <7q /3, the last asymptotic 
formula coincides with that given in (]55l) . 

Formula (1561) is also valid in the limit / — > go- However, 
in order to find explicit expressions for r and t 2 , it is much 
more convenient to use the basic formula (|36[) with / = go- 



A straightforward integration in this case yields 



vl 2(3-1 — 

go 2(p-iy T 



2(3- I 
2(/3 - 2) 



(59) 



for P > 1 and /3 > 2, respectively. Therefore, if / = go 
(i.e., g = l) and /3 > 2 then 



Dh 



igo P(J3 - 1) 



v (2/3-I) 2 (/3-2)' 



(60) 



A quite different behavior of Z?b on q occurs if the prob- 
ability density u(g) is mainly concentrated near the edges 
of the interval [— go, go] - In particular, for the limiting case 
characterized by the probability density 

u(g) = ^S(g-g ) + ^6(g + g ) (61) 



we obtain 



and so 



vl 



go <r 



l' 



Dh 



r2 - 



ho g 2 

2v a 



2 2 



W 2 



(62) 



(63) 



In contrast to the previous case, now Dh is a non-monoto- 
nic function of q (see figure [5]). According to (l6"3)l , Z?b ~ 
(lgo/2v)q~ 1 as g — > 00 (since in this case <r 2 = g 2 , the 
same result follows from (1551) as well), D^ ~ (lgo/v){q — 1) 
as g -> 1, and the maximum value of Z?b is -Db| 
0.192 x(lgo/v). 



7=V3 



5.3 Anomalous biased diffusion 

According to the above analysis, the biased diffusion is 
anomalous if both conditions / = go and a € (0, 2] hold. 
Specifically, if a £ (1,2] then r 2 = 00, r < 00 and so Z?b = 
00, i.e., the biased diffusion is fast in comparison with 
normal. In contrast, if a £ (0, 1] then r 2 = 00, r = 00, and 
depending on a the biased diffusion can be fast (Dh = 00), 
slow (Dh = 0) or even formally normal (0 < Dh < 00). 
Thus, taking into account that the points a = 1 and a = 2 
are distinct, we consider the long-time behavior of the 
variance a\(t) separately for a € (0, 1), a £ (1,2), a = 1, 
and a = 2. 



5.3.1 a £ (0,1). 

Since in this case the leading terms of (A t ) 2 and (Xf) are 
different (see below), for their finding we need to know 
only two leading terms of p s as s — > 0. These two terms 
can be easily evaluated using the expression 



Ps = 1 



dy(l 



>)P 



(64) 
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which follows from the definition p s — dre ST p(r) 
and the normalization condition drp(r) = 1 (we re- 
call that r m in = vl/2go and r max = oo at / = <?o)- Keeping 
in (|64p the leading term of the asymptotic expansion of the 
integral as s — > and using asymptotic formula (]37[) and 
the integral representation of the gamma function |33j . 
r(x) — J °° dy e~ y y x ~ 1 , we obtain the desired result 



Ps 



1 



ar(i 



(65) 



In order to find the leading terms of F s and G s , it 
is convenient to introduce the new variable of integration 
y = st and rewrite expressions (|18p and (]25[) in the form 



■y-i (v 



(1 + ^ + ^/2-1 { v_y 

y z \sJ 



and 



Using the asymptotic formula ([37]), one can make sure 
that the integral terms in and (1571) are proportional 
to l/s 1 "" and so F s ~ G s ~ 1/s as s — ► 0. Thus, in ac- 
cordance with this result and asymptotic expression f|65[) 
the leading terms of the Laplace transforms (IT9"]) and (|26p 
take the form 



la 



1 



a r(i - a) s x + a ' 



(xt 



2„,2 



2Pa 



1 



a 2 r 2 (l-a) s 1 + 2 "' 
(68) 

These asymptotic formulas are particular cases of the 
asymptotic formula (130[) in which the slowly varying func- 
tion L(l/s) is a constant. Therefore, as follows from (l3"Tj) , 
the long-time dependence of the moments (X t ) and (A' 2 ) 
is given by 



(Xt) 



Ic 



a,r(l-a)r(l + a) 



t° 



and 



(X 



(69) 



(70) 



a 2 r 2 (l -a)r(l+2a) 

As it seen from these asymptotic formulas, the leading 
terms of (A t ) 2 and (A 2 ) are different and so 



; 2 « 2 



a 2 r 2 (l - a) + 2a) r 2 (l + a) 



,2a 



(71) 

According to ([7T|) . the biased diffusion is characterized by 
a subdiffusive behavior (with = 0) if a G (0, 1/2) and 
by a superdiffusive one (with Z?b = oo) if a e (1/2,1). 
At a = 1/2 the diffusion is formally normal with = 
1 2 (tt — 2)/(27ra) 2 . However, since r = oo, r 2 = oo and 
(A t ) oc i 1 / 2 , we term this type of biased diffusion the 
quasi-normal diffusion. We note also that the asymptotic 
formula (1711) agrees with the asymptotic solution of the 
CTRWs OH]. 



5.3.2 a e (1,2) 



In this and all other cases the leading terms of (X t ) 2 
and (A 2 ) are the same. Therefore, for finding cr\{t) we 
need to determine three leading terms of p s as s — > 0. 
Since for these values of a the mean residence time f = 
f°°, dTTp(r) exists, to this end it is convenient to use the 
following exact formula: 



Ps 



1-ts + - s J^ dy(e- y +y-l)p(J) . (72) 
Proceeding in the same way as before, at s — > we obtain 

(73) 



_ ar(2-a) a 

Ps ~ 1 — TS H ; — S 



a(a — 1) 

Then, from (]5BJ and JBTJ it follows that F s ~ l/s-r/2 and 
Gs ~ 1/s— 2r/3 as s — > 0. Taking also into account that in 
the main approximation F s ~ 1/s and 2p s i ;l s + (l— p s )G s ~ 
2/s, the Laplace transforms (fTT)]) and (f2l))) are reduced to 



<*t>« 



Z 1 lar{2-a) 1 



^a(a - 1) s 3 



and 



2 2f_l_ 4:l 2 ar(2-a) 1 
( t)s ~ r 2 s 3 + r 3 a(a-l) s 4 - Q ' 



(74) 



(75) 



Finally, applying the modified Tauberian theorem to 
([74"]) and ([75]) and using the well-known property of the 
gamma function, _T(1 +x) = xr(x), we get 



(X t ) ~ = i + 



,2-a 



T 2 a{a - 1)(2 - a) 



and 



4Z 2 a 



»(a- l)(2-a)(3-a) 



(76) 



(77) 



as i — > oo. As it follows from ([76]) and ([77]). the long- 
time behavior of the variance a\ (t) obeys the asymptotic 
power law 



2Pa 



x3 — a 



r J a(2-a)(3-a) 



(78) 



It shows that at a € (1, 2) the unidirectional transport of 
particles is superdiffusive. The parameters in (|78j) depend 
on the probability density w(g) of the random force. In 
particular, if u(g) is determined by (]38[) then a = ft, the 
parameter a is given by (1391) . and 



ul a- 1/2 
5o a - 1 ' 



(79) 
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5.3.3 a = 1. 



5.3.4 a 



In order to illustrate the distinctive features of the long- 
time dependence of the moments (X t ) and (X 2 ) at a = 1, 
i.e., when according to ([57]) p(r) ~ a/r 2 (r — > oo), we first 
represent the Laplace transform p s in the form 



p s = 1 - aq s - 



dy(l - e~y) 



(80) 



where q s — s dy(l — e y )/y 2 - With the definition of 

the exponential integral [33], Ei(x) — dye~ v /y, the 
last quantity can be written as 



If a = 2 then it is convenient to represent the Laplace 
transform p s as 



1 — ts - 



where 



1 

ar s 

s 



dy(l-y-e v ) 



dy 



l-y-e 



1 - 2sr„ 



(1 



ST n 



2t£ 



(85) 



-Ei(sr min ). (86) 



1 - e 



Qs = 



sEi(sr min ). 



(81) 



Then, since p{r) — a/r 2 = o(l/r 2 ) as r — > oo (we recall 
that p(r) ~ a/r 2 if / = go and a = 1), the integral 
term in (|5U|) at s — > can be neglected in comparison 
with the term aq s . Finally, using the asymptotic formula 
£i(sT min ) ~ ln(l/s) (s -> 0) [33], we obtain 



p s 1 — as In — . 

s 



(82) 



As it can be easily seen, the integral terms in (|66D 
and ([67]) at a = 1 are proportional to ln(l/s) and so 
F s ~ G s ~ i/ s ( s — > 0). Therefore, the leading terms of 
the Laplace transforms (Ti~9")) and (l2"6l can be written in 
the form 



(Xt 



1 



a s 2 ln(l/s) ' 



(X 2 )s 



I 2 



a 2 s 3 ln 2 (l/s) 



(83) 



0) that in accordance with (|3Tj) yields 



(Xt 



I t 
a In t ' 



<* t 2 : 



a 2 hi 2 * 



(84) 



(t — > oo). Since (X t ) 2 ~ (^ 2 ), for finding the long-time 
behavior of the variance a 2 x (t) we need to know at least 
two terms of the asymptotic expansion of (Xt) and (X 2 ). 
However, in contrast to the previous case they cannot be 
determined from the modified Tauberian theorem because 
the leading terms of (X t ) s and (X 2 ) s contain the slowly 



varying functions, see (1831) . As it was mentioned in section 
15.2.11 in this case it is necessary to go beyond the Taube- 
rian theorem. Nevertheless, using the asymptotic laws (|71[) 
and (|78[) . we can make some qualitative conclusions about 
the character of biased diffusion at a = 1. First of all, it 
follows from and ((75} that t 2a and t 3 " a tend to i 2 
as a — > 1 — and a — > 1 + 0, respectively. But since 
r(l -a)->ooasa^l-0 and r— >ooasa— s-1 + 0, the 
coefficients of proportionality between cr^(t) and t 2 tend 
to zero. This means that at a = 1 the weakened ballistic 
diffusion should occur. 



Taking into account that p{r) — a/r 3 = o(l/r 3 ) (r — > 
oo), the integral term in ((85)) at s — >• can be neglected 
in comparison with ar s ~ — (as 2 /2) ln(l/s). Therefore, 
keeping in p s three leading terms, one obtains 



Ps 



1 



a 2, 1 
-s In-. 

2 s 



(87) 



Asymptotic expression (|87p shows that for determining 
two leading terms of (X t ) s and (X 2 ) s only the leading 
terms of F s and G s should be kept. Since F s ~ G s ~ 1/s, 
from ([TO]). (12^1) and ([57]) we immediately find 



lo 1, 1 
— — 7t — In-. 

2t 2 s s ' 



(X?) s 



2l 2 1 



2l 2 a 1 , 1 
^3 — 2 ln ~ 

T S z S 



(s —> 0). Accordingly, the modified Tauberian theorem, 
see (gni and (|1I]I. leads to 

I la I 2 2l 2 a 

(X t ) ~ -t+ -=olni, (AT 2 ) - =2 t 2 + ^-tlni (89) 
t 2r r T 

(i — > oo) and so the variance of the particle position in the 
long-time limit is described by the asymptotic formula 



(90) 



It shows that the biased diffusion at a — 2 is logarith- 
mically enhanced in comparison with normal diffusion oc- 
curring at a > 2. The fact that cr\(t) at t — > oo increases 
faster than t is in accordance with the diffusion law ([75]) . 
Indeed, while t 3 ~ a — > t as a —> 2 — 0, the coefficient of 
proportionality between the variance and time tends to 
infinity. 

For convenience, the above considered regimes of bi- 
ased diffusion that occurs under a constant force in a piece- 
wise linear random potential are summarized in figure [3] 
We note in this context that all laws of diffusion were ob- 
tained within a single model dealing with the solution X t 
of the motion equation @. In contrast, the CTRW jump 
model which deals with the discrete variable Y t reproduces 
only the long-time behavior of the variance erjf(t) if the 
waiting time is associated with the residence time [17] . At 
short times the behavior of a\ (€) and <j\ (€) is completely 
different, see sections E] 
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normal diffusion for all u(g) 



localization for all u(g) 







Fig. 3. Regimes of biased diffusion in a piecewise linear ran- 
dom potential. At / > go the normal biased diffusion occurs for 
all symmetric probability densities u(g) with support on the 
interval [—go, go]. A limiting case of subdiffusion, i.e., localiza- 
tion of particles at the average distance Q from the origin, 
is realized at / < go. If / = go then the character of diffu- 
sion depends on the exponent a characterizing the asymptotic 
behavior (|3Tf) of the residence time probability density p(r). 
Specifically, normal diffusion occurs at a > 2 (dashed line) 
and anomalous one at a € (0,2] (solid line). In turn, subdif- 
fusion occurs if a G (0, 1/2), superdiffusion if a £ (1/2, 1) or 
q G (1,2), and other types of anomalous diffusion are realized 
at a = 1/2 (quasi- normal diffusion), a — 1 (weakened bal- 
listic diffusion), and a — 2 (logarithmically enhanced normal 
diffusion). 



We complete our analysis of biased diffusion by com- 
paring the root-mean-square displacement ax(t) and the 
average displacement (X t ) at long times. Since (?x{t) char- 
acterizes the spreading of particles around their mean po- 
sition (Xt), the coefficient of variation C(t) = ax(t)/ (X t ) 
can be considered as a relative measure of the intensity 
of biased diffusion. Using the previous results, for nor- 
mal diffusion we obtain C(t) oc 1/t 1 / 2 (t — > oo). In the 
case of anomalous diffusion the coefficient of variation at 
a G [1,2] also tends to zero in the long-time limit, though 
more slowly than 1/i 1 / 2 . In particular, if a G (1,2) then, 
according to ([76) and ((78j), C{t) oc l/i (Q - 1)/2 . In contrast, 
as follows from (j6"9")l and (|7Tj) , the coefficient of variation 
C(oo) at a G (0, 1) does not vanish: 



C(oo) = 



2r 2 (l 



r{\ + 2a) 



- 1 



1/2 



(91) 



In this case the biased diffusion is so intensive that the 
spreading of particles is of the order of their displacement, 
i.e., ax(t) ~ C(oo)(X t ). Interestingly, since C(oo) — > as 
a — > 1 and C(oo) — > 1 as a — > 0, the slower the diffusion, 
the larger its intensity. 



Using the overdamped motion equation, we have derived 
the probability density of the particle position and have 
calculated its first two moments by the Laplace transform 
method. The Laplace transforms of these moments are 
represented by means of the probability density of the 
residence time, i.e., time that a particle spends moving on 
the interval of a fixed length, which in turn is explicitly 
expressed through the probability density of the random 
force. 

It has been shown that if the external force is less 
than the upper bound of the random force then the limit- 
ing case of subdiffusion, i.e., particle localization, occurs. 
In this regime, we have calculated the average distance 
between the origin and the points of localization which 
is always finite. In contrast, if the external force exceeds 
the upper bound of the random force then particles can 
be transported to an arbitrary large distance. Because of 
the influence of the random force, in this case the uni- 
directional motion of particles has a diffusive character. 
We have shown by applying the ordinary Tauberian the- 
orem for the Laplace transform that at small times the 
biased diffusion is always ballistic. For the analysis of the 
transport properties at long times we used the modified 
Tauberian theorem that in most cases permits us to find 
the two leading terms of the asymptotic expansion of the 
first and second moments of the particle position. Within 
this approach, we have shown that the biased diffusion is 
normal and the diffusion coefficient as a function of the 
external force can be either monotonic or non-monotonic. 
The latter occurs if the probability density of the random 
force is concentrated near the edges of the interval of sup- 
port. 

If the external force is equal to the boundary value 
of the random force then at short times the biased diffu- 
sion remains ballistic, but at long times it can be either 
normal or anomalous. In the last case the character of 
diffusion depends on the asymptotic behavior of the prob- 
ability density of the residence time, which is described by 
the exponent a. Using the modified Tauberian theorem, 
it has been shown that at a G (0, 2] the biased diffusion 
is anomalous. Specifically, subdiffusion, i.e., power-law de- 
pendence of the variance on time with power less than 1 
occurs at a G (0, 1/2), and superdiffusion, i.e., power-law 
dependence of the variance on time with power greater 
than 1 occurs at a G (1/2, 1) and a G (1,2). The quasi- 
normal diffusion, which is characterized by the nonlinear 
time dependence of the first moment of the particle posi- 
tion and the linear dependence of the variance, is realized 
at a = 1/2. Finally, the weakened ballistic diffusion and 
the logarithmically enhanced normal diffusion are realized 
at a = 1 and a = 2, respectively. 



6 Conclusions 

We have studied the unidirectional transport of particles 
which occurs under a constant external force in a piecewise 
linear random potential. The slopes of this potential, i.e., 
realizations of the piecewise constant random force, are 
assumed to be independent on different intervals of a fixed 
length and distributed with the same probability density. 
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